NaV 2 C>5 as a quarter-filled ladder compound 



Holger Smolinski, Claudius Gros, Werner Weber 
Institut fur Physik, Universitdt Dortmund, 44^21 Dortmund, Germany. 

Ulrich Peuchert, Georg Roth 
Institut fur Kristallographie, RWTH Aachen, 52056 Aachen, Germany 

Michael Weiden und Christoph Geibel 
MPI fur Chemische Physik fester Stoffe, TH-Darmstadt, 64289 Darmstadt, Germany 

(February 1, 2008) 



00 ■ 

o\ • 

ON i 

5h ; 

Oh 
<\ 

oo ; 

(N : 
> ■ 
\o ■ 

: 

(N . 

o ■ 

oo : 

ON . 

-i— > ■ 



C 

O 

o 



X 



A new X-ray diffraction study of the one- dimensional 
spin-Peierls compound a'-NaV^Os reveals a centrosymmetric 
(Pmmn) crystal structure with one type of V site, contrary to 
the previously postulated non-centrosymmetric P2imn struc- 
ture with two types of V sites (V +4 and V +5 ). Density func- 
tional calculations indicate that NaV2 0s is a quarter-filled 
ladder compound with the spins carried by V-O-V molecular 
orbitals on the rungs of the ladder. Estimates of the charge- 
transfer gap and the exchange coupling agree well with ex- 
periment and explain the insulating behavior of NaV2 0s and 
its magnetic properties. 

PACS numbers: 64.70.Kb, 71.27.+a, 75.10.-b, 75.30.Et 

Introduction Two classes of quasi one-dimensional 
compounds with spin-gaps have been investigated inten- 
sively in the last few years, ladder systems like SrCu203 
and Sri4Cu2404i [|]J|] and spin-Peierls compounds like 
CuGe0 3 and a'-NaV 2 5 H[§. For even-leg ladder sys- 
tems the spin-gap is structurally induced and present at 
all temperatures (for an overview see Q). The possibility 
of superconductivity in doped ladder systems has been 
discussed [[| and found in Sr2Cai2Cu2404i under pres- 
sure [p|. A spin-Peierls system undergoes a lattice insta- 
bility at Tsp and for T <T$p the system dimerizes and 
a spin-gap opens. Here we propose that the spin-Peierls 
compound NaV 2 Os is at the same time a quarter-filled 
ladder system, in contrast to the previous notion which 
assumed NaV205 to be made up of weakly coupled pairs 
of V +4 and V +5 chains [^j6) . Our proposition is based on 
a re-determination of the crystal structure of NaV20s by 
X-ray diffraction and on density- functional calculations. 
Mapping of the density functional results on Hubbard 
and Heisenberg models yields values for the model pa- 
rameters which explain readily the insulating behavior 
of NaV2 05 and its the magnetic properties. Our results 
also show that NaA^Os and CaV205 are isostructural 
and consequently establish CaV20s as a half- filled lad- 
der system. 

Crystal structure The crystal structure of a'-NaV^Os 
consists of double chains of edge-sharing distorted tetrag- 
onal VOs-pyramids running along the orthorhombic b- 
axis, which are linked together via common corners of 
the pyramids to form sheets. These in turn are stacked 



upon each other along c with no direct V-O-V links, (see 
Fig. [I] and j|j6j). The Na atoms are located between these 
sheets. In their original papers Galy et al. and Carpy 
et al. ^] proposed the non-centrosymmetric space group 
C2„-P2imn for this mixed valence compound (on average 
V +4 5 ). After the discovery of a spin-Peierls transition in 
a'-NaV 2 5 at T S p = 34 K f§, it has been argued that 
the charge ordering would lead to a magnetic decoupling 
of adjacent double chains j|,[7| and would be responsi- 
ble for the one-dimensional character of this compound 
observed in magnetic susceptibility measurements pLpfl. 
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FIG. 1. a) Crystal structure of NaV 2 5 . The dia- 
monds denote the position of the oxygen O(l)- and 0(2)-ions 
which form the basal quadrangles of the VO5 pyramids ori- 
ented in ±z directions. Also shown are the positions of the 
equivalent Vl/V2-atoms, located above/below the basal plan 
(empty/filled circles). Not indicated are the apex 0(3) and 
the Na positions. The dashed line indicates the intersection 
of the plane for the charge density plot of Fig. ||| b) Hopping 
matrix elements of our effective V(d xy ) Hamiltonian. 

Our re-determination of the structure by single crystal 
X-ray diffraction at room temperature, however, shows 
very clearly that the structure of NaV2 0s is in fact 
centrosymmetric (D^-Pmmn) with only one distinct V- 
position and 3 instead of 5 inequivalent oxygen atoms. 
NaV205 is consequently isostructural to CaV20s [p|,^0|. 
The relevant experimental details and results are summa- 
rized in Table |. The topology of the structure remains 
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essentially unchanged with respect to previous results || . 
The possibility for long range charge ordering, however, 
is lost in the higher symmetry group. 
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FIG. 2. V-d bands and density of states (DOS) as ob- 
tained from DFT. Also shown is a term scheme of the V-d 
orbitals, derived from an analysis of partial density of states. 
Here a,b,c,d, e denote respectively the d xy , d yz , d zx , d x i_ y i 
and d 3z 2_ r 2 orbitals. The filled dots in the second- lowest sub- 
band pair denote the range in k-space of 0(l)-p H admixture 
in the anti-bonding d xy band-pair. 

Attempts to refine the structure in the acentric space 
group P2imn resulted in unphysical anisotropic displace- 
ment parameters and strong correlations between posi- 
tional as well as displacement parameters of atoms which 
are equivalent in the centrosymmetric space group. De- 
spite the larger number of variables in the acentric space 
group, the reliability values did not improve at all. At- 
tempts to refine under the assumption of merohedral 
twinning (incoherent superposition of " inversion- twins" ) 
did not succeed either. Finally, a thorough inspection 
of a large number of centrosymmetric pairs of reflections 
(so-called Friedel-equivalents) did not reveal any signif- 
icant differences which could be indicative of a missing 
center of symmetry. These results are based on a number 
of measurements on various flux-grown crystals flH| from 
different batches, at two temperatures (293 K, 173 K), 
and with different wavelengths (MoKa, AgKa). Also, it 
was verified by Raman and susceptibility measurements 
p| that the crystals used for the diffraction experiments 
indeed exhibit a sharp spin-Peierls transition at 34 K. 
The data presented in Table || are for the homogeneous 
phase (T > T$p), the crystal structure for T <T$p has 
yet to be determined in detail Q] . 

Band structure Based on the new crystal structure 
(Table |) we have calculated the energy bands of NaV2 0s 
within density functional theory (DFT) jl^] . Thereby we 
have employed the full-potential linearized augmented 



plane wave code WIEN97 [|||. We have treated the 
exchange-correlation part by using the generalized gra- 
dient approximation |Q. Also, local orbitals have been 
included for a better description of the semicore states 
(of Na-2s, Na-2p, V-3s, V-3p and O-ls). 

The V-3d energy bands span a width of w 5 eV (see 
Fig. ||). The bottom of the 3d bands is separated by 
« 3 eV from the top of the 0-2p band manifold. The 
sequence of 3d subband-splittings is in accordance with 
estimates by ligand field theory from recent 51 V NMR 



data 1 15 



The four lowest-lying d bands predominantly exhibit 
V-d xy character (see Fig. |J). Actually, the d xy orbital 
planes are somewhat tilted around the b-axis towards 
the respective 0(1) positions above or below the centers 
of the V-O-V rungs (see Fig. ||). The four d xy bands are 
split into two pairs of subbands, separated by w 0.5 — lcV. 
All bands exhibit significant dispersion along T-Y, to a 
lesser extend along T-X, but hardly any dispersion along 
T-Z (see PigjJ)- 

Analysis jl6|j of the DFT band states and mapping of 
the bands on those of tight-binding models yield the fol- 
lowing result (model I): Bonding-type molecular orbital 
states, made up by the d^-orbitals of a V-O-V rung, 
coupled via the hopping term t± f=a —0.38 eV, form the 
lower pairs of subbands. Their dispersion along T-Y is 
produced by in « —0.17 eV, their splitting at T, the small 
dispersion along T-X and the band-crossing along T-Y re- 
sult from the small inter-ladder hoppings t\ ss 0.012 eV 
and *2 ~ 0.03 eV (terms see Fig. [l] b)). The upper pair 
of subbands consists of the corresponding anti-bonding 
molecular orbitals. 

In order to understand the microscopic origin of the 
effectiv hopping parameter t±, fn, t\ and t^., we have 
extended the tight-binding model by including the p x 
and p y orbitals of the basal plane 0(1) and 0(2) atoms 
(model II). We find strong contributions to t± and t\\ (of 
model I) by the indirect coupling of the V-d^y orbitals 
via the O-p, predominantly of pd7r-type. Further contri- 
butions to ij_ and in arise from residual three-center d-d 
terms involving the anionic oxygen potentials. The ratio 
ij_/i|| ~ 3.1 is related to the tilt of the V-d^ orbitals, 
which reduces the V-d xy — 0(2)-p x coupling along b 
and increases the Y-d xy — 0(l)-py coupling. The resid- 
ual direct d-d couplings are affected correspondingly by 
the tilt. 

Further, the indirect contribution to t\ via the VI- 
0(2)-V2 path is small due to a Goodenough-Kanamori- 
Anderson type interference effect in the (almost square- 
like) quadrangels of Vl-0(2)-V2-0(2), which is similar 
to that discussed in CuGeOa |l7]]. The residual direct 
contributions to t\ are affected by a compensation effect 
between two- and three-center contributions of opposite 
signs. 

The sign of in is derived from a symmetry analy- 
sis of the DFT eigenfunctions at the T-point. Admix- 
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ture of 0(l)-p y orbitals occurs only in the antibonding 
states of the V-O-V rung (upper subband pair). The p y - 
admixture is strongest near T, but disappears completely 
near Y (see Fig. ^). The latter feature points to a strong 
0(1)-0(1) pptr hopping (w 1 eV) along b. 
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FIG. 3. Partial charge density (e/A 3 ) of the V-d bands in 
the plane indicated Fig. This plane contains the Na sites 
and approaches closely the VI, V2 and apex 0(3) sites. The 
planes of the tilted d IH -orbital point 7r-like towards the O(l). 

The model resulting from our tight-binding analysis is 
that of a system of weakly coupled V-ladders with 0.5 
electrons per site (quarter filled). The band states are 
superpositions of d xy - d xy molecular orbitals of bonding 
and antibonding type. The pair of antibonding bands is 
empty and the the pair of bonding bands is half-filled. If 
we ignore the small inter-ladder couplings t% and t 2 , the 
d xy bond orbitals, which show up as the lowest subband- 
pair in Fig. ^|, form a pair of half-filled one-dimensional 
bands. Let us note, that this picture is robust and does 
not depend on the details of the tight-binding analysis. 
Inclusion of other intra-laddcr hopping parameters and 
0(l)-0(2) hopping matrix elements leads only to renor- 
malized tight-binding parameters, but not to a funda- 
mental modification of the model discussed above. De- 
tails will be discussed elsewhere. 

Correlation effects will lead to a charge gap and to 
the experimentally observed insulating behavior. Only 
one gapless spin-excitation branch remains |jl9|| . These 
statements hold both in the cases of weak and of strong 
electron-electron interaction U. We have estimated the 
value of U by a DFT-calculation where we have doped a 
fractional number of extra electrons. Though the doping 
is compensated on the Na sites, all additional electrons 
enter the V-d xy bands. The extra charge density causes 
a shift of all V-d bands with respect to the O-p bands. In 
Hartree-Fock approximation for our model I with on-site 
interaction U,the shift AE of the one-particle energies 
caused by the extra charge density An is given by jOl 

AE = f/An/2. 



For An = 0.1 we find a energy shift AE — 0.14 eV. We 
estimate U ~ 2.8 eV and a strong-coupling picture is 
therefore appropriate for NaV20s. We may thus regard 
NaV2 05 as being built up from weakly coupled spin-1/2 
chains with each spin located in a bonding V-O-V molec- 
ular wavefunction. Below T$p these chains of V-O-V 
spins dimerize and a spin-gap opens. 

t± — h~ 

o u— . . o 

'■ Jl 

-*±—t— 

FIG. 4. Term scheme of an empty, single and dou- 
bly-occupied V-O-V orbital representing an isolated rung of 
the ladder in the U S> t± limit. Given are the respective 
energies, J± » At\/U. 

We now proceed to estimate the charge-transfer gap 
E c and the antiferromagnetic coupling Jii. We neglect 
the inter-ladder couplings t\ and t 2 - Since t± ~ 3ty we 
may use perturbation-theory in tii /t±. The one- and two- 
particle states of an isolated V-O-V rung are depicted in 
Fig. |^. The bonding and antibonding one-particle states 
have the energies — t± and +t± respectively. The ex- 
change integral along a rung is J± = \ t\jU « 0.41 eV. 
In the ground state of the ladder all bonding states are 
filled and the first excited charge-transfer state is given 
by one empty and one doubly occupied rung with energy 

E c = 2tx - Ji ~ 0.71eV. 

This value for the charge-transfer agrees well with the 
0.6 — 0.7 eV observed in optical absorption spectra 
pot - The exchange coupling Ju between adjacent V-O-V 
molecular spins can be estimated by standard perturba- 
tion theory in t\\/t± and is given by 

2t 2 

J|| = -=r ~ 80 meV, 

which corresponds to 930 K. The exchange integral for 
NaV 2 5 has been estimated to be 560-700 K |Q. DFT 
therefore overestimates Jy somewhat. 

Outlook The above results indicate that ct'-NaV205 
may be the first known quarter-filled ladder compound 
It is possible to dope charge carriers into the in- 
sulating quarter-filled state. One way is to introduce Na 
defects, ^^205, as the a'-phase is stable for 0.7 < x < 
1.0 |22]. Alternatively one may consider the Ca substitu- 
tion, Nai_ y Caj / V 2 05, since our results establish NaV^Os 
and CaV2 05 to be isostructural |p|,|l0[| . Varying y e [0, 1] 
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would then allow to increase the carrier concentration 
continuously until a (highly anisotropic) half-filled ladder 
compound is obtained for y = 1. Note, that a spin-gap 
has been measured for CaN^Os |ic|]23| ]. 

Conclusions So far a'-NaX^Os has been considered 
as an inorganic spin-Peierls compound j||7|, where V +4 
and V +5 ions are ordered in parallel chains ||. In this 
letter we present evidence that NaV2C>5 is in fact a 
quarter-filled ladder system, consisting of equivalent V 
atoms, and that NaA^Os is isostructural to CaV205. 
Our results establish CaA^Os to be a half-filled lad- 
der compound and thus explain the observed spin-gap 
in CaV205. Our arguments are based on the crystallo- 
graphic re-examination of NaX^Os by X-ray diffraction 
and energy-band calculations. We find the crystal struc- 
ture to be centrosymmetric with only one equivalent V 
ion. This result does not allow for spontaneous charge 
disproportionation V +4 - V +5 . A tight-binding analysis 
of the band structure leads to a one-dimensional Heisen- 
berg model in the low-energy sector with a spin of 1/2 
per rung of the ladder. These spins are not attached to 
a single V ion, but to a V-O-V molecular orbital. Our 
estimates for the charge-transfer gap and the exchange 
coupling agree with experiment. 

We acknowledge useful discussions with P. Lemmens, 
R. Noack, A. Kluemper and M. Braden. 

Note added in proof: A recent V-NMR study also finds 
only one equivalent V-site in NaV20s above Tsp [ pi[ , as 
does a crystal structure redetermination The pos- 

sibility of a molecular spin-state has also been discussed 
recently p6l. 
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TABLE I. Data collection parameters: NaV^Os, Pmmn, 
compare Fig. |l|: a=11.316(4)A, b=3.611(l)A, c=4.797(2)A, 
z=2, T=293(2)K, A=0.56lA, sin6»/Aj max . = 1.365A" 1 , full 
sphere, numerical absorption correction, Ri n t= 0.052, 2279 
indep. data, 27 parameters, Rf=0.0238, wR|=0.0483 for 

I > 2(71. 

Atomic coordinates (xlO 4 ) and equivalent isotropic displace- 
ment parameters (A 2 x 10 3 ), see Table above. The two ad- 
ditional 0(2) depicted in Fig. [j] at (0.4269,0.75,0.5123) and 
(0.4269,-0.25,0.5123) are generated by symmetry. 
Bond lengths (A) : V-O(l): 1.8263(6), V-0(2): 1.9161(5) 
(2x), V-0(2): 1.9887(8) (lx), V-0(3): 1.6144(9), Na-O: 
2.5353(9) (avg.), shortest V-V: 3.0400(7) , shortest V-Na: 
3.3537(8). 
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